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ABSTRACT 

The relative abundances of ices in astrophysical environments rely on accurate laboratory measurements of physical 
parameters, such as band strengths (or absorption intensities), determined for the molecules of interest in relevant 
mixtures. In an extension of our previous study on pure-ice samples, here we focus on the near-infrared absorption 
features of molecules in mixtures with the dominant components of interstellar and planetary ices, H 2 0 and N 2 . We 
present experimentally measured near-infrared spectral information (peak positions, widths, and band strengths) for 
both H 2 0- and N 2 -dominated mixtures of CO (carbon monoxide), C0 2 (carbon dioxide), CH 4 (methane), and NH 3 
(ammonia). Band strengths were determined during sample deposition by correlating the growth of near-infrared 
features (10,000-4000 cm' 1 , 1-2.5 /Am) with better-known mid-infrared features (4000-400 cm' 1 , 2.5-25 /im) at 
longer wavelengths. 
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1. INTRODUCTION 

A variety of frozen volatile species (“ices”) are observed in the 
mantles coating the dust grains of the dense interstellar medium 
(ISM), where they are dominated by H 2 0 or N 2 (e.g., Gibb et al. 
2004), and on the surfaces of planetary bodies in the outer solar 
system (e.g., Roush 2001). Since the spectral properties (peak 
positions, spectral profile structures, and absorption strengths) 
of most species may be significantly altered by dilution, it is 
important to study them in an appropriate ice matrix. Here, we 
extend our previous study of the near-infrared spectral properties 
of pure ices (Gerakines et al. 2005) to those of molecules present 
within H 2 0- and N 2 -dominated ice mixtures. 

The compositions of interstellar icy grain mantles are dom- 
inated by H 2 0, by about a factor of two over the next-most 
abundant ice species, with minor components dependent on the 
specific environment probed along a chosen line of sight. Some 
of the most significant minor components of interstellar ices 
include CO (1%-10% relative to H 2 0), C0 2 (10%-30%), CH 4 
(l%-5%), and NH 3 (~1%; with some dependence on the envi- 
ronment sampled in the observed line of sight, e.g., Gibb et al. 
2004). Moreover, infrared spectroscopy of the vibrational ab- 
sorption features of these minor components indicates that while 
H 2 0-poor ice components do exist, a significant majority of icy 
mantle species appear to be diluted within an H 2 O-ice matrix 
(e.g., Tielens et al. 1991; Chiar et al. 1995; Boogert et al. 1996; 
Gerakines et al. 1999; Oliveira et al. 2009). 

In the outer solar system, many cold planetary surfaces appear 
or are thought to be dominated by H 2 0 (see, e.g., Hudson & 
Moore 2001), such as those of cometary nuclei (Roush 2001; 
Gibb et al. 2003), some Galilean satellites (e.g., McCord et al. 
1998; Hibbitts et al. 2000; Hansen & McCord 2008), many 
of the Saturnian satellites (e.g., Filacchione et al. 2007; Lopes 
et al. 2008), and many Kuiper Belt Objects (KBOs; Jewitt & Luu 
2004; Trujillo et al. 2005). Spectral modeling of the reflectance 
spectra indicate trace species on these surfaces such as CH 4 , 
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C0 2 , and NH 3 (e.g., Hibbitts et al. 2003; Cruikshank 2007; 
Schaller & Brown 2007; Clark et al. 2008). Molecules diluted 
in N 2 have been identified on planetary bodies such as Pluto and 
Triton (e.g., Cruikshank et al. 1998; Quirico et al. 1999; Buie & 
Grundy 2000; Grundy et al. 2003) and N 2 and CH 4 have been 
detected on large KBOs such as Sedna (KUO) (Barucci et al. 
2005) and Eris (Brown 2008). 

By and large, the current state of knowledge about the com- 
position of interstellar ices comes from observations and lab- 
oratory simulations of spectra obtained at mid-IR wavelengths 
(a = 2.5-25 /im). This is due to the fact that the strongest ab- 
sorption features of most small molecules in the solid phase fall 
in the mid-IR spectral region. These features are usually caused 
by fundamental molecular vibrational modes and are therefore 
the most sensitive tracers of ice composition (e.g., the 3.05 /Am 
O-H stretching vibration feature of H 2 0 or the 4.26 /Am C=0 
stretching vibration feature of C0 2 ). 

In contrast, the majority of what is known about solar system 
ices has been obtained from the spectra of reflected sunlight in 
the near-IR spectral region, where the absorption features due to 
overtones and combinations of the molecular vibrational modes 
are found. These features are weaker than the fundamental 
modes (e.g., Gerakines et al. 2005), but there is a higher flux 
of solar radiation in the near-IR. Most work in literature uses 
the wavelengths of absorption features in pure, unprocessed ices 
in comparison to observed data in order to deduce the general 
composition of icy surfaces (e.g., Roush 2001); however, recent 
experiments use ice mixtures (e.g., H 2 0 +NH 3 by Moore et al. 
2007 and N 1 +CH 4 by Brunetto et al. 2008). There is a need for 
laboratory mixtures in order to make accurate comparisons due 
to the fact that mathematical addition of individual spectra does 
not take into account the interactions of ice components when 
present in an intimate mixture. 

Studies of interstellar ices have extended into the near-IR 
spectral region as well. An example includes observations made 
with SpeX on the Infrared Telescope Facility telescope (Vacca 
et al. 2004). The James Webb Space Telescope Near Infrared 
Spectrometer will perform observations from 0.6 to 5 /Am at 
spectral resolutions from 100 to 3000 (Rauscher et al. 2004), 
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and will call for an increase in near-lR data. With our data, 
calculations of ice abundances in star-forming regions could be 
accomplished. 

A vital aspect of understanding ices on cryogenic surfaces 
in the ISM is the accurate determination of the abundances 
of the ice components. This relies on accurate laboratory 
measurements of physical parameters called band strengths. For 
example, the column density of a molecule in the solid phase ( N , 
in units of cm -2 ) as determined from astronomical observations 
of transmitted starlight is commonly calculated from the optical 
depth spectrum r(v) using the formula 


N= fi T (^dv 

At 


( 1 ) 


where the A, is the band strength of molecular absorption feature 
i (in units of cm molecule' 1 ), and the integration is taken over 
the wavenumber range of feature i. 

From Equation (1), it is clear that the accuracy of the compo- 
sitional information that may be derived from the astronomical 
data t(v) is closely tied to the accuracy of the band strength 
A,. Additionally, one cannot, in general, apply the same band 
strength to all ice mixtures containing the same molecular ab- 
sorption feature, but it must be determined for each mixture 
individually as this parameter is dependent on the specific phys- 
ical conditions of the local molecular environment (neighbor- 
ing molecules, crystalline structure, temperature, etc.; see, e.g., 
Gerakines et al. 1995 or Quirico & Schmitt 1997). The determi- 
nation of surface ice composition from a reflectance measure- 
ment also depends on the strengths of ice absorption features 
(or more precisely, the optical constants of the ice), as these 
parameters are required as input to surface scattering models 
used for interpretation. 

The IR absorption features of most molecules of astro- 
physical interest have band strengths from about 1 0 2 * * * " to 
10' 16 cm molecule” 1 , where the strongest absorption features 
generally lie in the mid-IR. Our laboratory group has already 
published a study of the near-IR band strengths of pure CO, 
CO 2 , C 3 O 2 , FhO, N'H ,, CH 4 , and CH 3 OH ice films at low tem- 
perature (Gerakines et al. 2005). This is a continuation of that 
work, using the same techniques, to the near-IR band strengths 
of selected mixed ices. Following the methods of Gerakines 
et al. (2005), here we present the near-IR band strengths of CO, 
CO 2 , CH 4 , and NH 3 as diluted in either N 2 or FFO, in the mixing 
ratio of ~5:1 at a temperature of 5 K. 


2. EXPERIMENTAL METHODS 

The experimental setup and ice creation techniques are iden- 
tical to those described by Gerakines et al. (2005) and simi- 
lar to those in previous studies by, e.g., Hudgins et al. (1993) 

or Gerakines et al. (1995). To summarize, the sample gases 

were deposited onto an IR-transparent substrate (Csl) at low 
temperature (~5 K) inside a high-vacuum sample chamber, 
and the resulting films were analyzed by FTIR spectroscopy. 
The vacuum system pressure was approximately 3 x 10 6 torr 
at room temperature before cooling the substrate. The sub- 

strate is attached to the end of a rotatable cold finger that 
may be cooled by a closed-cycle helium refrigerator (Air 
Products Displex 202N system). The temperature is continu- 
ously monitored by a chromel-Au thermocouple and is ad- 
justable by a resistive heater element from 5 K up to room 
temperature. 


The gases to be studied are mixed in the desired mixing ratio 
in a gas bulb attached to a separate manifold then transferred to 
the high- vacuum system. The sample gases are released within 
the vacuum chamber at a position about 1 cm in front of the 
substrate. Gas flow into the system is controlled by a needle 
valve, and flow rates are adjusted by this valve to achieve 
him growth rates of about 1-5 /im hr -1 . Two of the external 
ports of the vacuum chamber hold IR transparent windows 
made of KBr, and the entire chamber is situated inside the 
sample compartment of a ThermoMattson Infinity Gold FTIR 
spectrometer such that the IR beam passes at normal incidence 
through the KBr external ports and the Csl substrate. 

The spectrometer optics is capable of automatic switching 
between the configurations of sources, detectors, and beam 
splitters without breaking the dry air purge, allowing for 
the measurement of both near-IR spectra over the range of 
10,000-3500 cm -1 (1.0-2.86 /im) and mid-IR spectra from 
4500 to 400 cm ' 1 (2.22-25.0 /im). All spectra in this study 
were taken at a resolution of 4 cm -1 , which is suitable due to 
previous studies indicating the FWHM of these features was 
larger than 4 cm ' 1 (Gerakines et al. 2005). Ice samples created 
by vacuum deposition at ~5 K are expected to be amorphous 
rather than crystalline, and this was supported by the absorption 
profiles in the resulting IR spectra. The gases used and their 
purities are as follows: N 2 (Matheson, 98+%), CO (Matheson, 
99.9%), C0 2 (Matheson, 99.8%), CH 4 (Matheson, 99+%), NH 3 
(Matheson, 99+%), and H 2 O (distilled by freeze-thaw cycles 
under vacuum). 

By assuming a linear baseline, the integrated absorbance of 
each feature (in units of cm -1 ) was measured and the errors in 
the feature areas were estimated by integrating across the same 
limits when no feature was apparent above the noise level in 
the infrared spectrum. For narrow features (Au <10 cm' 1 ), the 
errors were found to be about 0.006-0.01 cm' 1 . This varied in 
direct proportion to a feature’s width. 

For a given molecule, the relative strengths of two absorption 
features were determined from the relationship between the 
areas of these bands as the ice was grown. While one could 
use a single ice spectrum and determine the relative strengths 
of any two features by taking the ratio of their areas, since 
the near-IR features are much weaker than those in the mid- 
IR (~10-100 times smaller; see absorbance scales in Figures 1 
and 2 ), errors in their measured areas are relatively much higher. 
Systematic errors involved in measuring the areas of the near-IR 
features are significantly reduced by monitoring the relationship 
between two given absorption features over the course of a slow 
deposit. 

In the experiments presented in the next section, we created 
mixtures of CO, CO 2 , CH 4 , and NH 3 diluted in N 2 and diluted 
in H 2 O with approximate compositions of dilutant: sample of 
5:1. These mixtures were chosen in order to represent some 
of the most general cases in both the ISM or the outer solar 
system. Dilution in H 2 O is meant to represent the environment 
of interstellar icy grain mantles or icy moons, such as Europa or 
Enceladus (Brown 2006; Hansen & McCord 2008). Dilution 
in N 2 is meant to represent the environment of icy moons 
dominated by N 2 , such as Triton (Buie & Grundy 2000; 
Stansberry 2004). Mixture ratios were not meant as a literal 
representation of astrophysical bodies, rather, mixtures were 
meant to be dilute enough to simulate a N 2 - or H 20 -dominated 
environment. 

For each sample molecule, the selected mixture was prepared 
in the gas manifold as described above. After cooling the 
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Figure 1. Combined near- and mid-IR spectra (! 0,000 400cm . 1.0-25 (im) of ice films at T ~ 5K for the N 2 -dominated mixtures. The sinusoidal baseline 

behaviors (strongest in the N 2 + CO 2 and N 2 + CO samples) are due to interference of the IR beam within the samples. Due to minor contamination by H 2 O and CO 2 , 
some features are also visible (see the text for a full discussion). 
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Figure 2. Combined near- and mid-IR spectra (1 0,000-400 cm” 1 , 1.0-25 /xm) of ice films at T ~ 5K for the H20-dominated mixtures. The sinusoidal baseline 
behaviors are due to interference of the IR beam within the samples. Due to minor contamination by CO 2 , some features are also visible (see the text for a full 
discussion). 
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Table 1 

Mid-infrared Features and Band Strengths Used in Strength Determinations 


Molecule 

V 

(cm -1 ) 

A 

W 

A 

(cm molecule -1 ) 

Reference 

Pure CH4 

1306 

7.657 

7.0 

X 

10- 18 

1 


2815 

3.552 

2.8 

X 

10- 19 

2 

CH 4 in H 2 0+ CH 4 (1.03:1) 

1302 

7.680 

8.1 

X 

10' 18 

3 

CH 4 in Nt+ CH 4 (10:1) 

2820 

3.546 

3.4 

X 

10' 19 

2 

Pure C0 2 a 

2283 

6.639 

7.8 

X 

10- 17 

4 


3708 

2.697 

1.4 

X 

10- 18 

4 

C0 2 in JDO+ C0 2 (1.6:1) 

2279 

6.639 

7.6 

X 

10' 17 

4 

Pure CO a 

2092 

4.780 

1.5 

X 

10-Wb 

5 

CO in H 2 0+ CO (2.1:1) 

2092 

4.780 

1.2 

X 

10 19b 

4 

Pure NH 3 a 

1070 

9.346 

1.2 

X 

10- 17 

5 

NH 3 in H 2 0+NH 3 (11:1) 

1070 

9.346 

1.3 

X 

10' 17 

3 


Notes. 

a Used for Ns mixture. See the text for further discussion. 
b Although this is a feature of 13 CO, its band strength is expressed in units of 
cm per 12 CO molecule. 

References. (1) Boogert et al. 1997; (2) Brunetto et al. 2008; (3) Kerkhof et al. 
1999; (4) Gerakines et al. 1995; (5) d'Hendecourt & Allamandola 1986. 

substrate to ~5 K, the gases were slowly bled into the vacuum 
system and the ice film was allowed to grow continuously 
over the course of several hours to days. The film growth 
was monitored at regular intervals, taking both near- and mid- 
IR spectra throughout the growth period. At each stage, the 
integrated areas of the mid- and near-IR features were measured. 
When it was determined that most of the near-IR absorption 
features were of sufficient strength to establish a clear trend, the 
gas flow was stopped and the experiment was ended. For each 
mixture involving FRO, it was necessary to replenish the gas 
bulb during the growth in order to deposit a sufficiently thick 
sample. 

3. RESULTS 

The spectra of ices containing solid CO, CO 2 , CH 4 , and NH 3 
diluted in N 2 in a ~5:1 mixture ratio and of ices containing solid 
CO, CO 2 , CH 4 , and NH 3 diluted in FbO in a ~5: 1 mixture ratio 
were collected in the near-IR region from 10,000 to 3500 cm ' 1 
( 1 .0-2.9 ;im) and the mid-IR region from 4500 to 400 cm 1 
(2.2-25 jim) during the slow growth of films at ~5 K. Figure 1 
displays the spectra of the ^-dominated samples and Figure 2 
displays the spectra of the HiO-dominated samples. If possible, 
the band strengths of any near-IR absorption features observed in 
the spectra were determined using the method described below. 

For optically thin ice absorptions, depositing a certain number 
of molecules should cause a linear increase the areas under all 
absorption features, where the rate of increase in area versus 
column density is equal to that feature’s band strength; see 
Equation (1). Hence, for a given molecule, the trend in one 
particular feature’s integrated area versus that of another feature 
of the same molecule should initially be a straight line whose 
slope corresponds to the ratio of the band strengths of the chosen 
features. In our case, we have chosen to compare the integrated 
areas of the observed near-IR features of each of the molecules 
under study to the areas of their much stronger mid-IR features. 
The mid-IR features used for these comparisons and their 
previously determined band strengths are listed in Table 1. 
Figures 3-6 display the near-IR regions of the spectra where 
absorptions were observed as well as the plots of the near-IR 
feature areas versus mid-IR feature areas for the ^-dominated 


ice mixtures. Figures 7-10 display the corresponding spectra 
and band area plots for the H 20 -dominated ice mixtures. For 
peaks where there is an underlining feature of H 2 O (see Figure 8 , 
for example), we drew a linear baseline underneath the feature 
to remove the underlying absorption due to H 2 O. 

Once the slope of the linear trend in the near-IR area versus 
mid-IR area is determined, it may be multiplied by the known 
band strength of the mid-IR feature (from Table 1) to obtain 
the near-IR band strength. Mid-IR band strength data exist for 
mixtures similar to the H 20 -ice mixtures and the N 2 + CH 4 
mixture studied here (Brunetto et al. 2008; Kerkhof et al. 1999; 
Gerakines et al. 1995) and those results were used to determine 
our near-IR band strength values. For CO, CO 2 , and NH 3 in 
N 2 -dominated mixtures, no such studies have been reported in 
the mid-IR. However, the studies of the molecules in H 2 O have 
shown that the mid-IR band strengths change by less than 20% 
relative to their pure-ice values (e.g., Gerakines et al. 1995). 
Since N 2 is a non-polar molecule, it can be expected that its 
effects would be even less pronounced. This is consistent with 
data for CH 4 (Kerkhof et al. 1999; Brunetto et al. 2008). Hence, 
where data for N 2 -dominated mixtures are unavailable, we have 
used the mid-IR band strengths for pure ices, with the caveat that 
this may introduce an error of approximately 20 % (but this is 
within the typical error in our measurements; see results below). 
Observed feature parameters (positions, widths, and vibrational 
assignments), slopes obtained for the linear trends in the near-IR 
versus mid-IR areas, and calculated band strengths using this 
method are given in Table 2 for N 2 -dominated mixtures and in 
Table 3 for H 20 -dominated mixtures. The listed uncertainties 
in near-IR band strengths were obtained by multiplying the 
standard deviations in the best-fitting slopes by the mid-IR band 
strengths and do not include the systematic errors associated 
with the band strength used for scaling. 

For some of the molecules (CO, CO 2 , and CH 4 ), the scaling 
process was complicated by the fact that some mid-IR absorp- 
tions became saturated (with peak absorbance values greater 
than 1.0) early in the ice film growth. The curves of growth 
were highly nonlinear. Thus, weaker features needed to be used 
in order to produce linear trends. For the CO ice deposits the 
absorption feature of 13 CO at 2092 cm 1 (4.780 /im) was used 
to scale the 12 CO and 13 CO near-IR features. The strength of 
1.5 x 10 -19 cm per 12 CO molecule was used for the 2092 cm ' 1 
(4.780 jim) feature, taking the band strength of the 13 CO feature 
from Gerakines et al. (1995) and scaling by the terrestrial ratio 
of 12 C/ 13 C of 87. For N 2 + CO 2 (5:1), the combination mode at 
3711 cm ' 1 (2.695 /xm) was used. For the H 2 O + CO 2 (5:1) ice 
deposit, the 2280 cm 1 feature of 13 C 02 was used in the same 
manner. For the H 2 O + CH 4 (5:1) ice the mid-IR absorption fea- 
ture at 1306 cm -1 (7.657 jim) was used for scaling (strengths 
are listed in Table 1). 

In the N 2 + CH 4 ice, the feature due to the CH 4 , the 
i >2 + V 4 vibration mode at 2825 cm -1 (3.540 jim), was used, 
assuming a band strength of A = 2.8 x 10 -19 cm molecule -1 
(corresponding to the 2815 cm -1 feature in pure CH 4 ice as 
determined by Brunetto et al. 2008). Note, in Gerakines et al. 
2005, a scaling error caused each of the pure CH 4 results to be 
off by a factor of seven and was corrected for future work. 

A small level of contamination due to CO 2 (~2%, as 
determined from the size of the absorption) was also observed 
within the N 2 ice mixtures, as is seen in Figure 1. Additionally, 
a very small amount of H 2 O was seen in the N 2 ice mixture 
(<1%). The value is small enough to not be seen in the 
Figure 1, due to the short time scale required for ^-dominated 
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Figure 3. Top: selected near-IR absorbance spectra of an N 2 + CH 4 (5:1) ice mixture at T ~ 5 K for the most prominent features studied. Bottom: integrated areas (in 
cm -1 ) of CH 4 features vs. the area (in cm' 1 ) of the 2825 cm ' 1 (3.540 /iml feature during deposition at T ~ 5 K. The solid lines are linear fits to the solid symbols, 
and data points shown as empty symbols were omitted from the fits (see explanation in the text). 


mixtures versus the H^O-dominated mixtures, which required 
changing of gas bulbs and additional pumping time between 
scans (as described in the Experimental Methods Section). Near- 
IR features occurring at ~7300cm -1 and ~5200cm -1 were 
observed in many EEO-dominated spectra (e.g., in Figure 2). 
These features were determined to be caused by small amounts 
of CO 2 and N 2 contamination within the FFO used to create the 
gas mixtures. This was not expected to significantly affect the 
band strength results. 


4. COMPARISONS WITH PREVIOUS STUDIES 

In this section, we compare these results to those found in the 
study of pure ices by Gerakines et al. (2005) and to other works 
in the literature. 


4.1. N 2 + CH 4 

For CH 4 diluted in N 2 , no visible features were detected for 
the (14 + 2v 3 + 14) and 3v 3 vibrational modes near 8780 and 
9060 cm' 1 , however, this may have been due to noise in the 
spectrum. All other vibrational modes showed an increase in 
peak wavenumber relative to pure CH 4 (ranging from 10 cm ' 1 
to 23 cm' 1 , respectively). All of the features were broader than 
in pure CH 4 (by 2-1 1 cm -1 ), with the exception of the (14 + 14) 
mode, which saw no broadening or sharpening of the peak, and 
the [vi] + [14] + [14] + [14] and the (2v 3 + 214) modes, which 
saw sharpening of the peaks and a decrease in the FWHM 
by 4 cm -1 and 7 cm -1 , respectively. Band strengths for the 
(v 3 + 14), (V 2 + v 3 ), 2 v 3 , and (14 + 2v 3 ) features (at 4311, 4539, 
6000, and 7503 cm -1 , respectively) measured were lower than 
those in the pure ice by a factor of 3-13. The (14 + 14), Oh + 14), 
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Figure 4. Top: selected near-IR absorbance spectrum of an N 2 + CO 2 (5:1) ice mixture at ~5 K for the most prominent features studied. Bottom: integrated areas (in 
cm -1 ) of CO 2 features vs. the area (in cm -1 ) of the 3711 cm -1 (2.695 fx m) feature during deposition at ~5 K. The solid lines are linear fits to the solid symbols, and 
data points shown as empty symbols were omitted from the fits (see explanation in the text). 
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Figure 5. Top: selected near-IR absorbance spectrum of an N 2 + CO (5:1) ice mixture at ~5 K for the most prominent features studied. Bottom: integrated areas (in 
cm -1 ) of CO features vs. the area (in cm -1 ) of the 2092 cm -1 (4.780 fim) feature during deposition at ~5 K. The solid lines are linear fits to the solid symbols. 
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Figure 6. Top: selected near-IR absorbance spectra of an N 2 + NH 3 (5:1) ice mixture at ~5 K for the most prominent features studied. Bottom: integrated areas (in 
cm -1 ) of NH 3 features vs. the area (in cm -1 ) of the 1 146 cm -1 (8.726 /xm) feature during deposition at ~5 K. The solid lines are linear fit to the solid symbols, and 
data points shown as empty symbols were omitted from the fits (see explanation in the text). 


Wavelength [|_im] 
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Figure 7. Top: selected near-IR absorbance spectrum of an H 2 O + CH 4 (5:1) ice mixture at ~5 K for the prominent most features studied. Bottom: integrated areas 
(in cm -1 ) of CH 4 features vs. the area (in cm -1 ) of the 1306 cm -1 (7.657 /xm) feature during deposition at ~5 K. The solid lines are linear fits to the solid symbols. 
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Figure 8. Top: selected near-IR absorbance spectrum of an H 2 O + CO 2 (5:1) ice mixture at ~5 K for the most prominent features studied. Bottom: integrated areas 
(in cm -1 ) of CO 2 features vs. the area (in cm -1 ) of the 2280 cm -1 (4.386 // m } feature during deposition at ~5 K. The solid lines are linear fits to the solid symbols, 
and data points shown as empty symbols were omitted from the fits (see explanation in the text). 
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Figure 9. Top: selected near-IR absorbance spectrum of an H 2 O + CO (5:1) ice mixture at ~5 K for the most prominent features studied. Bottom: integrated areas (in 
cm -1 ) of CO features vs. the area (in cm -1 ) of the 2092 cm -1 (4.386 // m ) feature during deposition at ~5 K. The solid lines are linear fits to the solid symbols. 


(113 + 204 ), and (V 2 + V 3 + V 4 ) features (at 4125, 4222, 5576, and 
5811 cm -1 , respectively) saw little to no change in the band 
strengths in comparison to the pure ice (<10%). Due to the 
nonlinear behavior of the growth curve at short deposition times, 
the strength of the [vi] + [vi] + [ 04 ] + [ 04 ] feature at 8422 cm -1 
could not be determined with accuracy. 

Near-IR spectra of N 2 + CH 4 in (1:1) and (10:1) ratios were 
reported by Brunetto et al. (2008), and our band strengths 
are comparable. For all features, our 5:1 mixtures indicated 
slightly larger peak shifts relative to the pure ice of an additional 
3-14 cm -1 . 


4.2. N 2 + CO 2 

For CO 2 diluted in N 2 , all expected near-IR vibrational 
modes were observed and showed a positive in peak position 
relative to the pure CO 2 ranging from 3 cnT 1 to 11 cm -1 , 
with the exception of the ( 2 vi + 2 t >2 + V 3 ) mode and the 3^3 
mode, which showed a decreases in position of 1 and 2 cm -1 , 
respectively. A slight broadening relative to pure CO 2 occurred 
in most features, ranging from 1 to 2 cnT 1 , with the excep- 
tion of the (t>i + 2 t >2 + V 3 ) and (2vi + P 3 ) modes, which showed 
no broadening or sharpening, and the (4ip + v 2 ) feature, which 
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Wavelength [|um] 

1.5 1.55 2 2.1 2.2 



Wavenumber [cm' 1 ] 



Figure 10. Top: selected near-IR absorbance spectra of an H 2 O + NH 3 (5:1) ice mixture at ~5 K for the most prominent features studied. Bottom: integrated areas (in 
cm -1 ) of NH 3 features vs. the area (in cm -1 ) of the 1110 cm -1 (9.009 /im) feature during deposition at ~5 K. The solid lines are linear fits to the solid symbols, and 
data points shown as empty symbols were omitted from the fits (see explanation in the text). 


showed a sharpening by 6 cm - 1 . All measurable band strengths 
showed an increase by 1.3-2 times those of pure CO 2 . Only 
small broadening and shifting occurred due to the nature of 
the apolar ices in the mixture (CO 2 and other apolar ices are 
not typically self-interactive). The increase in band strength 
may have been due to the slight contamination mentioned in 
Section 3. 


4.3. N 2 + CO 

For CO diluted in N 2 , all expected near-IR vibrational 
modes were observed and showed little change in peak po- 
sition relative to pure CO (+1 cm -1 for all modes, except 
the 2i> l 3 CO, which showed no change in position). Broad- 
ening relative to pure CO occurred in each feature ranging 
between 1 and 5 cm -1 . The only measurable band strength 
was that of the 2v I2 CO vibrational mode. The band strength 
was very comparable to the pure-ice feature’s band strength 
( 1.1 ± 0 . 1 ) x 10 -19 cm molecule -1 versus ( 1.6 ± 0 . 1 ) x 10 -19 
cm molecule -1 for the pure ice). As was seen with CO 2 , CO 
saw only minimal broadening and shifts and almost no change 
in the band strength due to the nature of two apolar ices in 
mixture, indicating almost no change in the interaction of the 
molecules. 


4.4. N 2 + NH 3 

For NH 3 diluted in N 2 , all expected near-IR vibrational modes 
were observed, with large shifts in the peak positions relative to 
pure NH 3 . Decreases in position occurred in both the (vj +v 2 + V 4 ) 
and the (v\+v 2 ) modes (of 65 and 30 cm - 1 , respectively), while 
increases were seen in the peak positions of the 2vi and (vj + V4) 
features (of 83 and 49 cm -1 ). All features showed a sharpening 
of peaks relative to pure NH 3 with decreases of 55 cm -1 to 
97 cm -1 . There was a decrease in the band strengths for the 
(vi + V2) and (vi + V4) features by a factor of 2^1 in comparison 
to those of pure NH 3 and an increase in the strength of the 2v\ 
mode feature by a factor of two. When diluted with N 2 , the 
mixture has fewer NH 3 -NH 3 interactions, likely resulting in the 
sharpening of the peaks. 

4.5. H 2 0 + CH 4 

For CH 4 diluted in H 2 O, no visible absorption features were 
detected for the (v 2 + 2 v 4 ) vibrational mode near 41 15 cm -1 . All 
other vibrational modes showed an increase in peak wavenum- 
ber relative to pure CH 4 (ranging from 4 to 213 cm -1 , respec- 
tively), with the exception of the (y 2 + V3) and 3 V 3 modes at 
4528 and 9039 cm -1 (the peak position of the 3 v 3 mode de- 
creased by 21 cm -1 from the pure-ice value). Most of the CH 4 
absorption features in the FFO mixtures were broader than their 


9 


The Astrophysical Journal, 759:74 (13pp), 2012 November 1 


Richey & Gerakines 


Table 2 

Observed Near-IR Features and Band Strengths of N 2 -dominated Ice Mixtures at 5 K 


Peak Position 
v (cm -1 ) 

X (/im) 

FWHM, Av 
(cm -1 ) 

Vibration 

Mode 

Relative Growth 11 
Rate 

Strength, A 
(cm molecule -1 ) 

N 2 +CH 4 (5:1) 






4125 

2.425 

14 

V2 + 2 v 4 

(3.9 ±0.3) x 10 -2 

(1.3 ±0.1) x 10 -20 

4222 

2.374 

11 

V\ + V4 

(79 ±4) x 10 -2 

(2.7 ±0.1) x 10 -19 

4311 

2.320 

18 

V 3 + V4 

(230 ±4) x 10 -2 

(7.9 ±0.1) x 10 -19 

4539 

2.203 

15 

V 2 + V 3 

(27 ±0.4) x 10 -2 

(9.2 ±0.1) x 10 -2 ° 

5576 

1.793 

15 

l>3 + 2V 4 

(3.6 ±0.4) x 10 -2 

( 1.2 ± 0 . 1 ) x 10 -2 ° 

5811 

1.721 

15 

V 2 + V 3 + V 4 

(7.1 ±0.4) x 10 -2 

(2.4 ±0.1) x 10 -20 

6000 

1.667 

18 

2 v 3 

(11 ±0.4) x 10 -2 

(3.6 ±0.1) x 10 -2 ° 

7503 

1.333 

18 

V 2 + 2 v 3 

(20 ±0.4) x 10 -3 

( 6.8 ± 0 . 1 ) x 10 -21 

8422 

1.187 

26 

Ft] + [i’ll + [14] + [V 4 ] 



8611 

1.161 

35 

2 v 3 + 2 v 4 

(5.4 ±0.4) x 10 -2 

( 1.8 ± 0 . 1 ) x 10 -2 ° 

N 2 +CO 2 (5:1) 






4840 

2.066 

7 

4v2 + V 3 

(9.3 ±0.7) x 10 -3 

(1.3 ±0.1) x 10 -2 ° 

4975 

2.010 

7 

v\ + 2 v 2 + V 3 

(59 ±0.7) x 10 -3 

( 8.2 ± 0 . 1 ) x 10 -20 

5090 

1.965 

7 

2 vi + V 3 

(25 ±0.7) x 10 -3 

(3.5 ±0.1) x 10 -2 ° 

6218 

1.608 

8 

V\ + 4v2 + V 3 



6340 

1.577 

9 

2 vi + 2 v 2 + ^3 



6970 

1.435 

10 

3v 3 



9955 

1.005 

6 

4 v 3 + V2(?) 



N 2 +CO (5:1) 






4159 

2.404 

7 

2v 13 CO 



4253 

2.351 

6 

2v 12 CO 

(73 ±7) x 10 -2 

( 1.1 ± 0 . 1 ) x 10 -19 

6339 

1.578 

9 

3v 12 CO 



8505 

1.176 

9 

? 



N 2 +NH 3 (5:1) 






4444 

2.250 

22 

V\ + V 2 (?) 

(18 ± 0 . 8 ) x 10 -3 

( 2.2 ± 0 . 1 ) x 10 -19 

5042 

1.983 

13 

Vl + 14 (?) 

(33 ±0.8) x 10 -3 

(4.0 ±0.1) x 10 -19 

6034 

1.657 


Vl + V >2 + ^4 (?) 



6598 

1.516 

21 

2 vi (?) 

(7.2 ±0.2) x 10 -3 

( 8.6 ± 0 . 2 ) x 10 -2 ° 


Notes. Vibrational mode assignments, except those in NH 3 mixtures, from Calvani et al. (1992) and Quirico & Schmitt (1997); NH 3 assignments are 
uncertain, however, are similar to those made in Moore et al. (2007). 
a The slope of the linear trend in the near-IR vs. mid-IR band areas shown in Figures 3-6. 


corresponding pure-ice values, with increases in FWHM from 
1 to 19 cm -1 , with the exception of the [vi] + [vi] + [V 4 ] + [ 114 ], 
(2u 3 + 2v 4 ), and 3v 3 modes (at 8569, 8612, and 9039 cm -1 ), 
which all showed a sharpening of the features with respect to 
their pure-ice values by 6-35 cm -1 . Due to the nonlinear be- 
havior of the growth curves even at short deposition times, the 
strengths of only two near-IR features could be determined with 
accuracy, as listed in Table 3, and each were found to be weaker 
than in the pure ice. 

AHiOt CH 4 mixture was studied at 15 K in two dilution 
ratios(20:l and3:l)by Bernstein etal. (2006) from 1.5 to 25 /xm 
(6700 to 400 cm -1 ). Generally, our measured widths and band 
strengths agree with those reported for the (v 3 + V 4 ), (v 3 + V 4 ), 
(V 2 + v 3 ), (i >2 + v 3 + V 4 ), and 2v 3 modes. For all features, our 5:1 
mixtures displayed slightly higher peak shifts than Bernstein 
et al. (2006) relative to the pure ice of an additional 4-10 cm -1 . 

Galvez et al. (2009) also studied several mixtures of FFO + 
CH 4 at 14 K. We find that our results agree in widths and band 
strengths for their H 2 O + CH 4 ( 10 : 1 ) and (4:1) ice mixtures for 
the ( vi +V4) and ( v 3 + V4) modes. Peak positions were also similar, 
with only a differences of 5-8 cm -1 . 

4.6. H 2 0 + CO 2 

For CO 2 diluted in H 2 O, no visible features were detected 
for the ( 4 v >2 + v 3 ) and (4v 3 + v 2 ) vibrational modes near 4840 


and 9955 cm -1 . The (vi + 2v 2 + v 3 ) at 4965 cm -1 decreased its 
peak positions by 6 cm -1 , respectively, in comparison to pure 
CO 2 value, while the (2vi + v 3 ) mode at 5087 cm -1 saw a slight 
decrease in position by 1 cm - 1 . Larger shifts were seen for the 
(v 3 + 4v 2 + v 3 ) mode at 6051 cm -1 (decrease of 163 cm -1 ), 
the (2v>i + 2 v 2 + v 3 ) mode at 6325 cm -1 (decrease of 16 cm -1 ), 
and the 3v 3 mode at 6864 cm -1 (decrease of 108 cm -1 ). While 
broadening occurred in every peak, the most prominent peaks 
saw the least amount of increase in FWHM. The (vi + 2v 2 + v 3 ) 
and (2vi + v 3 ) modes at 4965 and 5086 cm -1 saw broadening 
of only 2 and 3 cm -1 , respectively, while the (v\ + 4v 2 + v 3 ), 
(2vi + 2 v 2 + v 3 ), an d 3v 3 modes at 6051, 6325, and 6864 cm -1 
saw much larger increases in FWHM (93, 71, and 171 cm -1 , 
respectively). The strengths of only two absorption features 
could be determined with confidence: the (y\ + 2v 2 + n 3 ) and 
(2vj + 2 v 2 + u 3 ) modes, which showed a strength of 2-4 times 
those of the pure CO 2 band strengths. 

H 2 O + CO 2 was also studied by Bernstein et al. (2005) 
at 15 K in two mixtures, 25:1 and 5:1, from 2 to 25 /xm 
(5000-400 cm -1 ). Several near-IR bands were detected, includ- 
ing the (vi + 2 v 2 + v 3 ) feature, which in the 5:1 mixture showed 
a lower band strength 3.4 x 10 -20 cm molecule -1 compared to 
our value of (2.3 ± 0.1) x 10 -19 cm molecule - 1 and a larger 
FWHM (15 cm -1 compared to 9 cm -1 ) than in our results. In 
the 25:1 mixture of Bernstein et al. (2005), the band strength 
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Table 3 

Observed Near-IR Features and Band Strengths in H 20 -dominated Ice Mixtures at 5 K 


Peak Position 
v (cm -1 ) 

A. (/xm) 

FWHM, Av 
(cm” 1 ) 

Vibration 

Mode 

Relative Growth 3 
Rate 

Strength, A 
(cm molecule -1 ) 

H 2 O+CH 4 (5:1) 

4209 

2.376 

17 

Vl + V 4 

(3.6 ±0.1) x 10 ” 2 

(2.9 ±0.1) x 10 ” 19 

4308 

2.321 

16 

V 3 + V 4 

(7.0 ±0.3) x 10 ” 2 

(5.7 ±0.2) x 10 ” 19 

4528 

2.208 

10 

V 2 + f 3 



5640 

1.773 

8 

V 3 + 2 v 4 



5804 

1.723 

23 

V2 + V3 + V4 



5995 

1.668 

11 

2 v 3 



7700 

1.299 

26 

v >2 + 2 v 3 



8569 

1.167 

22 

[VI ] + [Vl] + [v 4 ] + [v 4 ] 



8612 

1.161 

7 

2 v 3 + 2 v 4 



8962 

1.116 

36 

v >2 + 2 v 3 + V 4 



9039 

1.106 

22 

3r 3 



H 2 O+CO 2 (5:1) 

4965 

2.014 

9 

vi + 2 v 2 + V 3 

(3.1 ±0.1) x 10 " 3 

(2.3 ±0.1) x 10 " 19 

5086 

1.966 

10 

2 vi + V 3 

(10 ±0.5) x 10 ” 4 

(7.6 ±0.4) x 10” 2 ° 

6051 

1.653 

100 

Vi + 4V2 + V 3 



6325 

1.581 

78 

2 vi + 2 v 2 + V 3 



6864 

1.457 

179 

3v 3 



H 2 O+CO (5:1) 

4260 

2.347 

33 

2v 12 CO 

1.2 ±0.7 

(1.4 ±0.8) x 10 ” 19 

6339 

1.575 

5 

3v 12 CO 



H 2 0+NH 3 (5:1) 

4518 

2.213 

22 

Vl + v 2 (?) 

(44 ±0.8) x 10 ” 3 

(5.7 ±0.1) x 10 ” 19 

5015 

1.994 

13 

Vl + v 4 (?) 

(38 ±0.8) x 10 ” 3 

(4.9 ±0.1) x 10 ” 19 

6136 

1.630 


Vl + V2 + V4 (?) 



6597 

1.516 

21 

2 vi (?) 

(27 ±0.8) x 10 ” 4 

(3.9 ±0.1) x 10” 2 ° 


Notes. Vibrational mode assignments, except those in NH 3 mixtures, from Calvani et al. (1992) and Quirico & Schmitt (1997); NH 3 assignments are 
uncertain, however, are similar to those made in Moore et al. (2007). 
a The slope of the linear trend in the near-IR vs. mid-IR band areas shown in Figures 7-10. 


and FWHM of this feature was consistent with those in our 
5:1 mixture. In each case, we measured slightly different peak 
positions. The (2vj + v 3 ) mode was not detected by Bernstein 
et al. (2005) in their 25:1 mixture; however, the properties of 
this feature in the 5:1 mixture were similar to those measured 
here. Bernstein et al. (2005) also observed two additional peaks, 
the 2v 3 and 4v 2 + v 3 modes, that were not visible in our spectra. 
This is likely due to overlapping H2O features (see Figure 2 or 
Figure 8). Another reason for this may be that the 2v 3 mode 
is a forbidden overtone of C0 3 and is extremely sensitive to 
conditions such as concentration and temperature. 

4.7. H 2 0 + CO 

For CO diluted in FFO, only the 2v 12 CO and 3u l2 CO features 
were visible at 4260 and 6339 cm” 1 , respectively. Compared to 
a pure ice, each showed an increase in peak position (by 8 cm' 1 
and 1 cm -1 , respectively) and each broadened, by 30 cm -1 
and 1 cm -1 , respectively. Only the strength of the 2v I2 CO 
feature was able to be determined from our data ((1.4 ±0.1) x 
10” 19 cm molecule” 1 ), and it was found to be comparable to the 
strength in the pure ice ((1.6 ± 0.1) x 10” 19 cm molecule” 1 ). 

4.8. H 2 0 +NH 3 

For the NH 3 diluted in FFO, all observed features were 
seen to shift in peak position in comparison to pure NFI 3 , 
with increases ranging from 22 to 82 cm” 1 , and to decrease in 
FWHM by 55-97 cm” 1 . The strengths of all detectable features 
also decreased from their pure-ice values by approximately 


1.5 times. The strength of the (v 3 + 2v 2 + V4) mode near 
6136 cm” 1 could not be determined from our data. 

H 2 0 + NH 3 (—3:1) at 9 K was studied by Moore et al. (2007). 
Their results included peak positions of the (fi + v 2 ) and (ui + 114 ) 
vibrational modes. For the (v t + v 2 ) feature, Moore et al. (2007) 
reported a peak position of 4478 cm” 1 , which is further (by 
7 cm” 1 ) from the pure-ice value than our results indicate. The 
(v\ + V 4 ) mode also had a larger change from the pure-ice value 
than seen in our results (by 2 cm” 1 ). Moore et al. (2007) did not 
calculate bands strengths. 

5. DISCUSSION AND IMPLICATIONS 

These experimental data can be used for more powerful in- 
terpretations of observations of astrophysical ices. In particular, 
the band strengths measured here may be used to more accu- 
rately determine column densities of ice components in the ISM. 
As shown in Gerakines et al. (2005), using our laboratory data, 
we can predict the optical depths of near-IR ice absorptions 
for several astrophysical objects. We have estimated the optical 
depths for several lines of sight in the ISM using the widths 
and strengths (as well as previously published column densi- 
ties) for some of the near-IR features of our H20-dominated 
mixtures (e.g., Gerakines et al. 1999; Gibb et al. 2004). These 
estimates are listed in Table 4 and, due to their direct depen- 
dence on band strengths, show weakening in comparison to 
the results predicted in Gerakines et al. (2005). The estimates 
here are more appropriate to interstellar ices than our pure-ice 
analysis in Gerakines et al. (2005), since interstellar ices are 
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Table 4 

Column Densities 11 and Predicted Near-IR Optical Depths’ 3 of Interstellar Ices Diluted in H 2 O 



W33A 

NGC 7538 IRS9 

Elias 16 

Sgr A* 

GL 7009S 

GL 2136 

ch 4 

N 

r(2.376 /xm) 
r(2.321 /xm) 

1.7 x 10 17 
0.003 
0.006 

1.5 x 10 17 
0.002 
0.006 


2.6 x 10 16 
<0.001 
<0.001 

4.8 x 10 17 
0.008 
0.01 


co 2 

N 

r(2.014/rm) 
r(1.966 (im) 

1.4 x 10 18 
0.04 
0.01 

1.7 x 10 18 
0.05 
0.01 

4.5 x 10 17 
0.01 
0.004 

1.8 x 10 17 
0.005 
.001 

2.5 x 10 18 
0.07 
0.02 

8.0 x 10 17 
0.02 
0.006 

CO 

N 

r(2.347 fim) 

8.8 x 10 17 
0.004 

1.2 x 10 18 
0.005 

6.3 x 10 17 
0.002 

<2 x 10 17 
<0.001 

1.8 x 10 18 
0.008 

1.0 x 10 17 
<0.001 

nh 3 

N 

r(2.213 /xm) 
r(1.994 /xm) 

1.7 x 10 18 
0.04 
0.06 

9.8 x 10 17 
0.02 
0.03 

<2.3 x 10 17 
~ <0.006 
<0.009 

2.6 -3.9 x 10 17 
0.006-0.01 
0.01-0.02 




Notes. 

a Column densities (/V) have units of cm -2 ; Gibb et al. (2004). 

b Optical depths have been calculated from the listed column densities, using width and strength data listed in Table 3. 


dominated by water ice. However, they may not fully reflect 
the most realistic cases for these objects, as some are located 
in highly energetic environments, such as NGC 7538 IRS 9, 
GL 7009S, and GL 2136, which are all in the vicinity of young 
stellar objects and therefore undergoing various degrees of heat- 
ing and/or energetic processing (Gibb et al. 2000). 

N 2 -dominated ice spectra may be used in comparison with ob- 
served planetary ices instead of pure ices to deduce information 
about the general composition of the surface ice. An example is 
Figure 2 in Quirico et al. (1999) in which the spectrum of Triton 
is compared with laboratory transmission spectrum of pure ices 
CH 4 , CO, N 2 , CO 2 , and H 2 O. Similar to Triton, Pluto has been 
recently investigated by Buie et al. (2010) and has shown signs 
of N 2 ice, however, the behavior (i.e., seasonal effects) of ices 
on both planetary surfaces remains poorly understood. With the 
upcoming New Horizons mission in 2015, the data presented 
here may be useful for interpretation and further understanding 
of the general composition of the planetary surface. However, 
the chemical compositions of icy planetary bodies are affected 
by particle radiation or high-energy photons, resulting in the 
synthesis of new molecular species which can then be ejected 
into the gas phase by various mechanisms (Bringa & Johnson 
2004). Studying the effect of UV radiolysis on these ices in ad- 
dition to placing them in molecular mixture would lead to more 
realistic choices for comparison to observational data. 

6 . SUMMARY AND FUTURE WORK 

In this paper, we have collected spectra of N 2 mixed with 
CO, CO 2 , CH 4 , and NH 3 in 5:1 ratios and of H 2 O mixed 
with CO, CO 2 , CH 4 , and NH 3 in 5:1 ratios in the near-IR 
region from 10,000 to 3500 cm -1 (1.0-2. 9 /im) and the mid- 
IR region from 4500 to 400 cm ” 1 (2.2-25 /im) during the slow 
growth of films at ~5 K and have used these to the determine 
near-IR band strengths listed in Tables 2 and 3. Future work 
will analyze the effect of UV radiolysis on both the pure ices 
previously discussed in Gerakines et al. (2005) and the ice 
mixtures described in this paper. We will also investigate the 
effects, if any, of ice composition on the near-IR band strengths 


as well as the differences between the near-IR band strengths of 
crystalline and amorphous ices. 
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